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I review the theory of preheating after inflation, focussing on the recently 
VO \ found tachyonic preheating in the theories with spontaneous symmetry break- 

ing. This occurs due to the tachyonic instability of the scalar field near the 
top of its effective potential. Contrary to the common expectation, tachyonic 
instability converts most of the energy into that of colliding classical waves 
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very rapidly, within a single oscillation. Efficient tachyonic preheating is typi- 
cal for the hybrid inflationary scenario, including SUSY motivated and brane 
^^ | inflation models. 

O; 

^ ■ I. PREHEATING AFTER INFLATION 

a: 

Ph! The best-fit universe is the uniform and homogeneous flat Q tot = 1 Friedmann-Lemetre 

expanding universe with gaussian scale-free metric fluctuations. The model is in remarkable 
agreement with the observations of the CMB anisotropy, the large scale structure of the 



X 



universe and the tests of the global geometry. 

Inflationary theory suggests that the expansion of the universe was preceded not by the 
Big Bang singularity, but rather by the exponentially fast expansion of the universe when 
four dimensional geometry of the universe is close to the nonsingular De Sitter geometry. 
Almost De Sitter geometry (where the Hubble parameter is slow varying with time, H <C H 2 ) 
can be realized by gravity of the matter with the vacuum-like equation of state p ~ —p, 
where p and p are the pressure and energy density. There are many specific realizations of 
this equation of state, motivated by the different aspects of the elementary particle physics. 
The vacuum-like equation of state can be provided by the out-of-equilibrium scalar fields or 
by the other effects (like the quantum gravity mechanisms) which behave like an effective 
scalar field. 

Although scalar fields are not yet discovered experimentally, they are the vital ingredients 
of the high energy physics theories, and a plethora of scalar fields exists in the supergarvity 
and superstring theories. Fundamental M-theory should encompass both supergravity and 
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string theory. At present the low-energy phenomenology is described by the N = 1 d 
supergravity. Its Lagrangian begins with the scalar fields terms 
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where a scalar field $ l is the complex conjugate of $j. Some preferable choices of the Kahler 
potentials K, superpotentials W and Yang-Mills couplings hopefully will be selected at the 
level of the fundamental theory. Until the fundamental theory of all interactions is well 
understood, one may try to address the issues of the early universe cosmology in the context 
of the most general phenomenological N = 1 supergravity- Yang-Mills-matter theory. This, 
in fact, was the case during the last decade. 

At first glance it seems hard to gain useful insights or predictions from a loose pool of 
inflationary models embedded in (|l|). Fortunately, theory of inflation, independently on its 
concrete model, provides us with several universal mechanisms which make it immensely 
successful as the theory of the initial conditions for the observable Friedmann-Lemetre uni- 
verse. These mechanisms are related to the properties of the DeSitter solution in General 
Relativity and to the properties of classical and quantum fields in the background of De 
Sitter space-time. Remarkably, these properties are not directly related to the microscopic 
physics behind the inflation. This is not the first time when we can successfully do cosmol- 
ogy at the phenomenological level without the knowledge of the microscopic physics. For 
example, celebrated Cold Dark Matter theory of the Large Scale Structure of the universe 
relies on the simple assumption that there is a nonbaryonic invisible dark matter component 
with the "dust-like" equation of state p c dm ~ 0. This is enough to develop the theory of the 
structure formation even without understanding the microscopic origin of dark matter. 

There are several major predictions of inflation: 

• Inflation erases pre-inflationary classical inhomogeneities and entropy. 

• The total mass density in the universe Qtot = 1- 

• Vacuum quantum fluctuations of the scalar field(s) generate almost scale free gaussian 
scalar metric perturbations. 

• Vacuum quantum fluctuations of gravitational waves generate almost scale free gaussian 
tensor metric perturbations. 

• All the particles in the universe are created from the decay of inflaton energy p in the 
process of (p)reheating after inflation. 

The ratio of amplitudes of scalar and tensor modes, the spectral index of scalar and tensor 
modes, and the character of (p)reheating after inflation can be model dependent. 

In this contribution I will review the theory of preheating after inflation, concentrating 
on the recently found tachyonic preheating [|10| , |TT|1 . 

According to the inflationary scenario, the Universe initially expands quasi-exponentially 
in a vacuum-like (DeSitter) state without entropy or particles. At the stage of inflation, all 
energy is contained in a classical slowly moving fields $ in the inflaton sector. The last term 
of the 1st line of ([[]) is the scalar potential V($j). The equations of motion based on the first 
line should describe inflation, which is a challenging problem by itself. The Lagrangian ([!]) 
contains also other fields which give subdominant contributions to gravity. The Friedmann 
equation for the scale factor a(t) and the equation for <&(£) determine the evolution of the 
background fields. 



In the chaotic inflation models, soon after the end of inflation, an almost homogeneous 
inflaton field $(£) coherently oscillates with a very large amplitude of the order of the Planck 
mass Mp around the minimum of its potential V{&). This scalar field can be considered 
as a coherent superposition of inflatons with zero momenta. The amplitude of oscillations 
gradually decreases not only because of the expansion of the universe, but also because 
energy is transferred to particles created by the oscillating field. At this stage we shall recall 
the rest of the fundamental Lagrangian which includes all the fields interacting with inflaton. 
These interactions lead to the creation of many ultra- relativist ic particles from the inflaton. 
Gradually, the inflaton field decays and transfers all of its energy to the created particles. In 
this scenario all the matter constituting the universe is created from this process of reheating. 
If the creation of particles is sufficiently slow, the particles would simultaneously interact 
with each other and come to a state of thermal equilibrium at the reheating temperature 
Tr. Typically particle production from coherently oscillating inflatons occurs not in the 
non-perturbative regime of preheating |T]|. 

Indeed, let us consider a simple toy model of chaotic inflation with the quadratic potential 
V($) = |m^$ 2 and Ci n t = — \g 2 $ 2 X 2 describing the interaction between the inflatons 
$ = $1 and other massless Bose particles \ = $2- We can consider quantum fluctuations of 
the field x interacting with the classical homogeneous background field $(£). The quantum 
scalar field x in a flat FRW background has the eigenf unctions Xfc(^) e ~ ikx with comoving 
momentum k. The temporal part of the eigenfunction obeys the equation 
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Xk + 3-Xk + K - & + 9 2 4> 2 }Xk = (2) 
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with vacuum-like initial conditions: Xk — ^/sr in the far past. The coupling to the curvature 
£R will not be important in the presence of the interaction (but would lead to gravitational 
preheating in the absence of the interaction). In this model, the inflaton field $(£) coherently 
oscillates as $(£) ~ $(£) sin(m^t), with the amplitude $(£) = -j== • -^— t decreasing as the 
universe expands. 

Equation (||) describes excitation of the quantum fluctuation Xk- At first glance the effect 
of particle creation n x could be treated perturbatively with respect to the small coupling g 2 . 
However, the smallness of g 2 alone does not necessarily lead to the perturbative approach 
to describe the excitation of Xk modes. To check whether the interaction term g 2 <p 2 in eq. 
(0) is perturbative or not, we have to use a new time variable z = mt and the essential 
dimensionless coupling parameter q = a — 3-. Scalar metric fluctuations in this model are 
compatible with cosmology if the inflaton mass is m^ ~ 10~ 6 M p ; therefore, typically q ~ 
10 10 g 2 ^> 1 for not negligibly small g 2 . In fact, a consistent setting for the problem of 
X-particle creation from the 0-inflaton requires g> 1 even without additional assumptions 
about g 2 . It is known that if we have two scalars <j) and Xi then the latest stage of inflation 
will be driven by the lightest scalar. The square of the effective mass of the x-Aeld includes 
a term g 2 <f) 2 . Inflation is driven by the 0-field if its square mass m 2 is smaller than g 2 <f) 2 . 
This leads to the condition q ^> 1. 

Solutions of the equation (||) for q ^> 1 are unstable for a range of k and x particles are 
created exponentially fast. Indeed, suppose that there is no expansion of the universe. In 
this case background oscillations are harmonic, $(£) ~ sin (m^t), and equation (0) is reduced 
to the Mathieu equation. Solutions of this equation are exponentially unstable within the 



set of resonance bands, Xk ~ e^" 1 ^*, where the characteristic exponent/ifc depends on (k,q). 
In realistic case of the expanding universe, parameter q is time dependent. For the broad 
resonance case q ^> 1 this parameter jumps over a number of instability bands within a 
single background oscillation, so the concept of stability /instability bands is inapplicable 
here. Parametric resonance in this case, described by (0), is a stochastic process M. In 
the regime of stochastic resonance particle are created exponentially fast as n^ ~ e 2 J * Mfc . 
Characteristic exponent here is also stochastic with typical value flk ~ 0.1 in the resonance 
momentum range < k < q l ^m. 

Due to the copious particle creation from the background field, very soon, within few 
dozens of the background oscillations, we have to take into accound backreaction of the cre- 
ated particles. Unfortunately, known theoretical approaches to include backreaction effects, 
like the Hartree approximation, are not enough in this situation. First fully nonlinear lattice 
simulation of the preheating || revealed that the leading backreaction effect is rescattering 
of x particles on the background inlatons, X^o — > Xp^k-k'- Lattice simulations demon- 
strated that there is a first, resonant stage of preheating, where initial excitation of x field 
occurs exponentially fast as it is expected from the linear theory of equation (Q). Then the 
system enters the nonlinear stage of preheating where classical waves of x an d $ fields are 
rescattering and gradually relaxing towards equilibrium. 

Different in technical details, but qualitatively similar theory was developed for other 
types of the chaotic inflation potentials, e.g. see [Q for V($) = |A$ 4 theory. 

II. PREHEATING IN HYBRID INFLATION 

Another popular class of inflationary models - hybrid inflation - involve multiple scalar 
fields $j in the inflaton sector ||. In particular, hybrid inflation can be realized in su- 
pergravity for certain choice of superpoteantial W 0, and in the theory of brane inflation 

0- 

Previous studies of preheating in hybrid models were concentrated on particle creation by 
parametric resonance that may occur when homogeneous background fields oscillate around 
the minimum of the potential @,§]. Technically the linear fluctuations were considered 
around the time- dependent homogeneous background fields. Such parametric resonance may 
or may not be strong depending on the coupling parameters. However, we recently found 
Til that there is strong preheating in hybrid inflation, but its character is quite different 



from preheating based on parametric resonance. It turns out that there is very efficient 
tachyonic instability that appears in the hybrid inflation models. The backreaction of rapidly 
generated fluctuations does not allow homogeneous background oscillations to occur because 
all energy of the oscillating field is transferred to the energy of long-wavelength scalar field 
fluctuations within a single oscillation! However, this does not preclude the subsequent decay 
of the Higgs and inflaton inhomogeneities into other particles, and thus reheating without 
parametric resonance. 

Let us first consider the background evolution and the results of the naive perturbative 
approach to decsribe the quantum fluctuations around the background solutions. Consider 
the simple potential for the two-field hybrid inflation is 

V{^a) = h^-v*y + 9 ^o\ (3) 



where we used notations $1 = 0, $2 = o. Inflation in this model occurs while the homoge- 
neous $1 field slow rolls from large <fi towards the bifurcation point at <\> = —v (due to the 
slight lift of the potential in direction). Once 0(t) crosses the bifurcation point, the cur- 
vature of the a field, m 2 = d 2 V/da 2 , becomes negative. This negative curvature results in 
exponential growth of a fluctuations. Inflation then ends abruptly in a "waterfall" manner. 

One reason to be interested in hybrid inflation is that it can be implemented in super- 
symmetric theories ([[]). In particular, for illustration we will consider preheating in the 
super symmetric F-term inflation as an example of a hybrid model. 

The simplest F-term hybrid inflation model (without undesirable domaine walls) is based 
on a superpotential with three left-chiral superfields $j = ($1, $2, $3) in the Lagrangian ([[]) 

W = ^-$1 (4$ 2 $3 - v 2 ) . (4) 

In this case, the spontaneous breaking of the local (global) U(l) symmetry between the $2 
and $3 fields will lead to gauge (global) string formation. 

In global SUSY, using the same notation for superfields and their complex scalar com- 
ponents, this superpotential contributes 

V F = ^|4$2$3 - v 2 \ 2 + 4A|$!| 2 (|$ 3 | 2 + |$ 2 | 2 ) • (5) 

to the scalar potential. In general, $3 and $3 could be (oppositely) charged under a local 
£7(1) symmetry, in which case we should include a D-term, Vd, which we neglect here. 

In this model, inflation occurs when chaotic initial conditions lead to (|$i|) 3> v. When 
this happens, the fields $2 and $3 acquire large effective masses and roll to their local 
minimum at ($2) = (^3) = 0. In this limit, the potential (d) becomes V ~ ^-, which 
gives rise to a non- vanishing effective cosmological constant. However, this is a false vacuum 
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state; the true vacuum corresponds to ($2 ( I ) 3) = \ and ($1) = 0. The slow-roll potential 
drives the evolution of the inflaton towards its true VEV. When its magnitude reaches the 
value (|$ c |) — § spontaneous symmetry breaking occurs. 

For further discussion of symmetry breaking in this model, let us rewrite ([5]) in terms of 
polar fields: $3 = |$3|e* e , $2 = |$2|e* 9 . The potential becomes 

V F = -(l6\$ 2 \ 2 \& i \ 2 -8v 2 \$ 2 \\$ 3 \cos(8 + d)+v 4 ) + 4A|$ 2 | 2 (|$ 2 | 2 + |$ 3 | 2 ) • (6) 

At the stage of symmetry breaking, when (|$2$3|) begins to move away from zero, the 
absolute phase Arg($ 2 $3) = + 9 acquires a mass and is forced to zero. Note, however, 
that the potential is independent of the relative phase, 8 — 8, reflecting the U(l) symmetry. 
Thus, in a quasi- homogeneous patch, the U(l) symmetry allows us to choose the relative 
phase of the $2, ^3 fields to be zero without any loss of generality. This choice, combined 
with the vanishing of the absolute phase, is equivalent to choosing the two complex $ 2 , 
$3 fields to be real. In order to leave canonical kinetic terms, we define a± = |$3| ± l^l- 
Furthermore, as inflation has left the inflaton homogeneous across all the patches, we may 
choose it to be real: <p = \/2|$i|. 

In terms of these three real fields, the potential now becomes 



V F = j (oj - a 2 - v 2 f + A0 2 (a 2 + + a 2 ) . (7) 

In the symmetric phase, when a± = 0, the a fields have an effective mass-squared: m±((j)) = 
A (2(p 2 =F v 2 ). We can now see that spontaneous symmetry breaking occurs in this model 
exactly as in the two field model (|3]). For (f> < <f> c = -t=, the <r + field has a tachyonic mass 
that triggers symmetry breaking and the end of inflation. On the other hand, the <r_ field has 
always a large and positive effective mass-squared, pinning it to zero. Thus, during inflation 
and at the initial stages of symmetry breaking, this model behaves just like the standard 
two field hybrid model (^). We have only to apply the constraint g 2 = 2A and identify the 
Higgs field with a + . The equations for the homogeneous background components <p(t) and 
u±(t) admits simple solution 

0(t) + -^x + (£)=0 c , <r_(t) = 0. (8) 

To study preheating in the F-term inflation, we have to analyse evolution of the vacuum 
fluctuations. Consider vacuum fluctuations in the inflaton sector $j of the theory (H). Usual 
description in terms of a homogeneous background plus small fluctuations gives us equations 
for fluctuations around the background solution (||). We define the variances of fields: 

°± ~ (°±) I 2 ) rcn = / (03 [l&M*)| 2 - l^ ± (0)| 2 ] =Jj-V ± {k 7 t) (9) 

and similar for field. Here V±^(k 1 t) are the spectra of the fluctuations. 
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FIG. 1. Evolution of the background fields o+(£) and <f)(t) after symmetry breaking and the log 
of V+(k) for the mode with the momentum k = 0.2vAf • 
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FIG. 2. Means and variances in units of <p c . The squared means ((ft) 2 and (<r+) 2 are ordinary 
solid lines while the field variances are thick lines. The mean of (f> starts at <j) c , oscillates once, and 
then decays. The mean of 0"+ grows in antiphase to <p and freeses at <p c . 

Numerical solutions of the equations for the background fields and for the time evolution 
of the linear fluctuations of a+ for the mode k where their spectrum is maximum are plotted 
as the bold line at the Figure [|. Notice an enormous exponential growth of the fluctuations 
within a single background oscillation. Indeed, the amplitude V+(h) increases by factor 
10 10 ! There are two factors which contributes to such a strong instability of fluctuations 
in the model. First, oscillating background fields are crossing the region with significant 
negative curvature of the effective potential, which results in tachyonic instability. Second, 
this region turn out to be a turning point for the background oscillations, where the fields 
spend significant portion of the oscillation. As a result, tachyonic instability is lasting long 
enough to make the backreaction of the fluctuations to be significant already within single 
background oscillation. The regime of background oscillations will even not be settled. 
Therefore practically from the beginning we have to use the lattice simulations to study 
nonlinear dynamics of the fields. 

The results of full nonlinear lattice simulations in the model derived in [jTUJ] are plotted 
in Figure 0. The simulations showed that the homogeneous fields <fi and a + initially followed 
the classical trajectory (§) but, within one oscillation of the inflaton field, fluctuations grew 
too large to speak meaningfully of the fields as homogeneous oscillators. These fluctuations 
grew in such a way that <j+ = cr_* almost exactly throughout the simulation. In other words 



Re 5a + and Im 5<j- were excited while Im 5a + and Re 5 a _ were not. Because of this we 
only plot the fields and a + . 

In fllt)| , |nf we develop a general theory of tachyonic preheating, which occurs due to 
tachyonic instability in the theories with spontaneous symmetry breaking. Our approach 
combines analytical estimates with lattice simulations taking into account all backreaction 
effects. The process of spontaneous symmetry breaking involves transfer of the potential 
energy into the energy of fluctuations produced due to the tachyonic instability. We show 
that this process is extremely efficient and requires just a single oscillation of the scalar field 
falling from the top of the effective potential. In what follows I will illustrate tachyonic 
preheating using the model (f|). 

III. SPONTANEOUS SYMMETRY BREAKING AND TACHYONIC 

PREHEATING 

To understand physics of tachyonic preheating, we have to go back to the basics of 
the spontaneous symmetry breaking [pj| . The simplest model of spontaneous symmetry 



breaking is based on the theory with the effective potential Vp = |f 4 + \a\ — ka\ + Aa 2 ., 
dominated by the quadratic term near the top of the potential. The development of tachyonic 
instability in this model depends on the initial conditions. We will assume that initially the 
symmetry is completely restored. Initially scalar field fluctuations in this model in the 
symmetric phase = are the same as for a massless field, <fik = -L^e~ lht+tkx . Then at 
t = we 'turn on' the term —m 2 <f) 2 /2 corresponding to the negative mass squared — m 2 . The 
modes with k — \k\ < m grow exponentially. Initial dispersion of all growing fluctuations 

with k < m was given by (5(j) 2 } = J |^- = |^ and the average initial amplitude of all 

fluctuations with k < m was given by 5(f) = ^ The dispersion of the growing modes at 
t > is growing exponentially. This means that the average amplitude 5<p(k) of quantum 
fluctuations with momenta ~ k initially was 5<f>(k) ~ k/2ir, and then it started growing 

as e tv ^ m ~ k . The tachyonic growth of all fluctuations with k < m continues until J (5(f) 2 ) 

reaches the value ~ v/2, since at ~ i>/v3 the curvature of the effective potential vanishes 
and instead of tachyonic growth one has the usual oscillations of all the modes. This happens 
within a time At ~ tt- In t-, where C ~ 10 2 . 

The process of symmetry breaking will occur in a somewhat different way in theories 
where the curvature of the effective potential near its maximum depends on 0. Consider for 
example V = — ^0 4 near the top of the potential. In this case there is an instanton solution 
corresponding to the tunneling from = to nonzero value [|3~2| , p~3|] . 

Let us now consider the tachyonic instability for the theory (|j). The scalar fields potential 
(0) for the solution (|J) in the fields space is reduced to the potential 

^ = -/ - \*l + ~a4 (io) 

Let us simplify notation a + — > and make rescaling A — > A/3. Then the potential ( |10|) once 
again is reduced to the form 

V = -^ + ^ + ±v\ (11) 



Thus the theory ([11]) is a prototype of the theory ([|), which initially looked much more 
complicated. 

Development of instability in the theory (|II]) presents us with a new challenge. The 
curvature of the effective potential at = in this theory vanishes, which means that, 
unlike in the theory — m 2 2 , infmitesimally small perturbations in this theory do not grow. 
On the other hand, in this theory, unlike in the theory — A0 4 , there are no instantons which 
would describe tunneling from = 0. Thus, in the theory —\v<p 3 , which occupies an 
intermediate position between — m 2 2 and — A0 4 , both mechanisms which could lead to the 
development of instability do not work. Does it mean that the state = in this theory is, 
in fact, stable? 

The answer to this question is no, the state = in the theory — Xv 3 is unstable. 
Indeed, even though (0) initially is zero, long wavelength fluctuations of the field are 
present, and they may play the same role as the homogeneous field in triggering the 
instability. 

The scalar field fluctuations with momentum k < k have initial amplitude (50 2 ) ~ 

k 2 

g^-. Thus the short wavelength fluctuations with momenta k > ko live on top of the long 
wavelength field with an average amplitude <50 rms (/co) ~ J (50 2 ) ~ oZk~- 

The curvature of the effective potential V" = |m 2 ff | at ~ S(p rms (k ) in the theory ( |Tl|) is 
given by — 2XvS(p rras (ko) ~ — Xv-M-. Consider the fluctuations with momentum k somewhat 
greater than ko, so that the amplitude of the long wavelength field 50 does not change 
significantly on a scale k~ x . Short wavelength fluctuations with k = Cko with C somewhat 
greater than 1 will grow on top of the field ~ <50 r ms(^o) if k 2 < |m 2 ff | ~ ^ a . 

Taking for definiteness C ~ v2, one may argue that fluctuations with k < j 1 may enter 
a self- sustained regime of tachyonic growth. Small fluctuations rapidly grow large, which 
justifies semi-classical methods used for the description of this process. The average initial 
amplitude of the growing tachyonic fluctuations with momenta smaller than ^ is 
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(12) 



These fluctuations grow until the amplitude of 50 becomes comparable to 2u/3, and the 
effective tachyonic mass vanishes. At that moment the field can be represented as a collection 



of waves with dispersion J(5(j) 2 } ~ v, corresponding to coherent states of scalar particles 
with occupation numbers n^ ~ f 4 ^ - ) ^> 1- A more accurate investigation shows that the 



initial value of the field is few times greater than J(5(j) 2 ) ~ v, and therefore the occupation 
numbers will be somewhat smaller, 

n fc ~O(10)A" 2 . (13) 

Because of the nonlinear dependence of the tachyonic mass on 0, a detailed description 
of this process is more involved than in the quadratic theory. Indeed, even though the 
typical amplitude of the growing fluctuations is given by (|T2"D, the speed of the growth 
of the fluctuations increases considerably if the initial amplitude is somewhat bigger than 
(|I2|). Thus even though the fluctuations with amplitude a few times greater than ( |T2"D are 
exponentially suppressed, they grow faster and may therefore have greater impact on the 
process than the fluctuations with amplitude (|12|). 






FIG. 3. Three panels above show growth of quantum fluctuations of the field 4> i n the theory 
V = —-^vcf) 3 + 4 4 (|ll| ) looks like bubble formation. Preheating occurs due to a combined effect of 
bubble production, tachyonic instability and bubble wall collisions. 



10 



Low probability fluctuations with £0 3> £0rms correspond to peaks of the initial Gaussian 

distribution of the fluctuations of the field 0. Such peaks tend to be spherically symmetric 

TM. As a result, the whole process looks not like a uniform growth of all modes, but more 



like bubble production (even though there are no instantons in this model). 

To study this issue in a more detailed way, one may use a stochastic approach to tunneling 
and bubble formation developed in |H| . The main idea of this approach can be explained as 
follows. Tunneling can be represented as a result of accumulation of quantum fluctuations 
with the amplitude greatly exceeding their usual value determined by uncertainty principle. 
This happens when the long wavelength quantum fluctuations responsible for the tunneling 
correspond to bosonic excitations with large occupation numbers. In such cases one can treat 
these fluctuations as classical fields experiencing Brownian motion due to their interaction 
with the short wavelength quantum fluctuations. 

Suppose that the large fluctuations of the scalar field responsible for reheating in the 
model ( P]) initially look like spherically symmetric bubbles (which is the case if the prob- 
ability of such fluctuations is strongly suppressed, see above). Equation of motion for a 
bubble of a scalar field 0(r) in Minkowski space is 



= 0" + 20V" 1 - V'{<j>) . (14) 

Here r is a distance from the center of the bubble, 0' = -S. At the moment of its formation, 
the bubble wall does not move, = 0, = (critical bubble). Then it gradually starts 
growing, > 0, which requires that 

|0" + 20V- 1 ! < -V'{<S>) . (15) 

A bubble of a classical field is formed only if it contains a sufficiently large field 0, and if the 
bubble itself is sufficiently large. If the size of the bubble is too small, the gradient terms 
are greater than the term |V"(0)|, and the field inside the bubble does not grow. 

At small r the shape of the bubble can be approximated by = 0(0) — ar 2 /2. In this 

approximation, the bubble has a typical size r ~ J -^, and 0'r -1 = 0" = —a. Therefore 
at the moment of the bubble formation, when = 0, one has 

0" = H0(O))/3 • (16) 

Replacing 0" by &q0(O) one finds that the bubble can be considered a result of overlapping 
of quantum fluctuations with typical momenta k < ko ~ r^ , where 

*J^^#. (17) 

30(0) V ; 

Here C = 0(1) is some numerical factor reflecting uncertainty in our estimate of k . 

Let us estimate the probability of an event when vacuum fluctuations occasionally build 
up a configuration of the field satisfying this condition. In order to do it one should remember 
that the dispersion of quantum fluctuations of the field with k < k is given by (Sep 2 ) ~ ^. 
This gives 

U 2, K 2 n0(o)) (m , 

{(P )k<k0 ~8^ = C 24^0(0) ■ (18) 
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This is an estimate of the dispersion of perturbations which may sum up to produce 
a bubble of the field that satisfies the condition (|15|) . Of course, this estimate is rather 
crude. But let us nevertheless use eq. ([18|) to evaluate the probability that these fluctuations 
build up a bubble of a radius r > k^ 1 containing the field at its center. Assuming, in the 
first approximation, that the probability distribution is gaussian, one finds: 

m ~ exp hfc) = exp {-z£Q • (l9) 

The general formula ([19]), being applied to the theory — A0 4 /4 to within a factor of C ~ 2 
coincides with the Euclidean action for the instanton in this theory. Taking into account the 
very rough method we used to estimate k and calculate the dispersion of the perturbations 
responsible for tunneling, the coincidence is rather impressive. As it was shown in [|16| , |15|1 , in 
application to the tunneling during inflation in the potentials with V" <C H 2 this approach 
gives exactly the same answer as the Euclidean approach, ost importantly, this methods 
allows to investigate tunneling and development of instability in the theories where the 
instanton solutions do not exist |I5| . In particular, for the tunneling in the theory — At>0 3 /3 
one finds 

P(0)~(A^) 2 exp(-^) . (20) 

We included here the subexponential factor 0(/c 4 ) ~ (Av0) 2 , which is necessary to describe 
the probability of tunneling per unit time per unit volume. 

This means that the tunneling is not suppressed for ~ j^ ■ This result is in agreement 



with our previous estimate ([12]) . Now let us take into account that the total time of the 
development of instability is a sum of the time of tunneling plus the time necessary for rolling 
of the field down. One can show that the time of rolling down is inversely proportional to 
m(0) ~ y/\vcj), i.e. it decreases at large 0. Also, the subexponential factor (At>0) 2 grows 
at large 0, which makes tunneling to large faster. Consequently, as we already discussed 
above, the main contribution to the development of instability is given by the fluctuations 
with >~ %^r- Exponential suppression of the probability of such fluctuations leads to 
their approximate spherical symmetry. 



The results of our lattice simulations for this model |l Ij are shown in three panel of Fig. 
|3[. In this model bubbles form quickly enough that we were able to start with quantum 
fluctuations centered at < >= and allow the bubbles to form automatically. The 
bubbles (high peaks of the field distribution) grow, change shape, and interact with each 
other, rapidly dissipating the vacuum energy V(0). 

IV. DEVELOPMENT OF EQUILIBRIUM AFTER PREHEATING 

The character of preheating may vary from model to model, e.g. parametric excitation 
in chaotic inflation |l| and tachyonic preheating in hybrid inflation [n|, but its distinct 



feature remains the same: rapid amplification of one or more bosonic fields to exponentially 
large occupation numbers. This amplification is eventually shut down by backreaction of 
the produced fluctuations. The end result of the process is a turbulent medium of coupled, 
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inhomogeneous, classical waves far from equilibrium. Despite the development of our un- 
derstanding of preheating after inflation, the transition from this stage to a hot Friedmann 
universe in thermal equilibrium has remained relatively poorly understood. The details of 
this thermalization stage depend on the constituents of the fundamental Lagrangian ([!]) and 
their couplings, so at first glance it would seem that a description of this process would 
have to be strongly model-dependent. Recently we performed a fully nonlinear study of 
the development of equilibrium after preheating |)17|| . We have performed lattice simula- 
tions of the evolution of interacting scalar fields during and after preheating for a variety 
of inflationary models. We have found, however, that many features of this stage seem to 
hold generically across a wide spectrum of models. Indeed, at the end of preheating and 
beginning of the turbulent stage £*, the fields are out of equilibrium. We have examined 
many models and found that at £* there is not much trace of the linear stage of preheating 
and conditions at £* are not qualitatively sensitive to the details of inflation. We therefore 
expect that this second, highly nonlinear, turbulent stage of preheating may exhibit some 
universal, model-independent features. Although a realistic model would include one or 
more Higgs- Yang-Mills sectors, we treat the simpler case of interacting scalars. 

We have numerically investigated the processes of preheating and thermalization in a 
variety of models and determined a set of rules that seem to hold generically. These rules 
can be formulated as follows (in this section we use notations <fi = $1 for the inflaton field 
and X) c for other scalars $;) 

1. In many, if not all viable models of inflation there exists a mechanism for exponen- 
tially amplifying fluctuations of at least one field x- These mechanisms tend to excite long- 
wavelength excitations, giving rise to a highly infrared spectrum. 

The mechanism of parametric resonance in single-field models of inflation has been stud- 
ied for a number of years. This effect is quite robust. Adding additional fields (e.g. a fields) 
or self-couplings (e.g. x 4 ) nas little or no effect on the resonant period. Moreover, in many 
hybrid models a similar effect occurs due to tachyonic instability. The qualitative features 
of the fields arising from these processes seem to be largely independent of the details of 
inflation or the mechanisms used to produce the fields. 

2. Exciting one field \ is sufficient to rapidly drag all other light fields with which \ interacts 
into a similarly excited state. 

We have seen this effect when multiple fields are coupled directly to x an d when chains 
of fields are coupled indirectly to \. All it takes is one field being excited to rapidly amplify 
an entire sector of interacting fields. These second generation amplified fields will inherit 
the basic features of the x field, i.e. they will have spectra with more energy in the infrared 
than would be expected for a thermal distribution. 

3. The excited fields will be grouped into subsets with identical characteristics (spectra, 
occupation numbers, effective temperatures) depending on the coupling strengths. 

We have seen this effect in a variety of models. For example in the models ( pip which 
we are going to consider the \ and a fields formed such a group. In general, fields that 
are interacting in a group such as this will thermalize much more quickly than other fields, 
presumably because they have more potential to interact and scatter particles into high 
momentum states. 
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4- Once the fields are amplified, they will approach thermal equilibrium by scattering energy 
into higher momentum modes. 

This process of thermalization involves a slow redistribution of the particle occupation 
number as low momentum particles are scattered and combined into higher momentum 
modes. The result of this scattering is to decrease the tilt of the infrared portion of the 
spectrum and increase the ultraviolet cutoff of the spectrum. Within each field group the 
evolution proceeds identically for all fields, but different groups can thermalize at very 
different rates. 

Here we will illustrate these results with a simple chaotic inflation model with a quartic 
inflaton potential. The inflaton <fi has a four-legs coupling to another scalar field x-, which 
in turn can couple to two other scalars o~\ and o-l. The potential for this model is 

V = \\<fi* + ^W + \h\ X 2 o\ + \h\ X 2 ol (21) 

_ 2 

Preheating in this theory in the absence of the <Tj fields is well studied. For nonsmall 3 r- 
the field x wm experience parametric amplification, rapidly rising to exponentially large 
occupation numbers. In the absence of the x field (or for sufficiently small g) <fi will be 
resonantly amplified through its own self-interaction, but this self-amplification is much 
less efficient than the two-field interaction. The results shown here are for A = 9 x 10 -14 
(for CMB anisotropy normalization) and g 2 = 200A. When we add a third field we use 
h\ = 100g 2 and when we add a fourth field we use h\ = 200g 2 . 

One of the most interesting variable to calculate is the (comoving) number density of 
particles of the fields n(t) and their occupation number n k . The evolution of the total 
number density of all particles n to t is an indication of the physical processes taking place. 
In the weak interaction limit the scattering of classical waves via the interaction \g 2 (j) 2 x 2 
can be treated using a perturbation expansion with respect to g 2 . The leading four- legs 
diagrams for this interaction corresponds to a two-particle collision ((fix ~^ 0x)> which con- 
serves ntot- The regime where such interactions dominate corresponds to "weak turbulence" 
in the terminology of the theory of wave turbulence. If we see n to t conserved it will be an 
indication that these two-particle collisions constitute the dominant interaction. Conversely, 
violation of n tot (t) = const will indicate the presence of strong turbulence, i.e. the impor- 
tance of many-particle collisions. Such higher order interactions may be significant despite 
the smallness of the coupling parameter g 2 (and others) because of the large occupation 
numbers n k . Later, when these occupation numbers are reduced by rescattering, the two- 
particle collision should become dominant and n to t should be conserved. For a bosonic field 
in thermal equilibrium with a temperature T and a chemical potential fi the spectrum of 
occupation numbers in the limit of classical waves is given by 

n k « . (22) 

L) k - fl 
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FIG. 4. Time evolution of number density of particles in the model (|2l|). The curves represent 
n^, n x , n (J1 , n ff2 from top to bottom. Unit of (conformal) time is a • 10 -36 sec. 

Figure [| shows an exponential increase of n{t) during preheating, followed by a gradual 
decrease that asymptotically slows down. This exponential increase is a consequence of 
explosive particle production due to parametric resonance. After preheating the fields enter 
a turbulent regime. In our simulations we see n(t) decreasing during this stage. This decrease 
is a consequence of the many-particle interactions beyond the four legs rescattering. 

An important point is that the interaction of x an d &i does not affect the preheating 
of Xii but does drag <Tj exponentially quickly into an excited state. The fields <jj are expo- 
nentially amplified not by parametric resonance, but by their stimulated interactions with 
the amplified \ field. Unlike amplification by preheating, this direct decay nearly conserves 
particle number, with the result that n x decreases as <Tj grow. 

Interacting waves of scalar fields constitute a dynamical system. Dynamical chaos is one 



of the features of wave turbulence. In |L7[] we address the question how and when the onset 
of chaos takes place after preheating. To investigate the onset of chaos we have to follow 
the time evolution of two initially nearby points in the phase space. Consider the theory 
with the potential (|21~|) with two fields <fi and \ only (which we collectively denote as /). 
Consider two configurations of a scalar field / and /' that are identical except for a small 
difference of the fields at a set of points xa- Chaos can be defined as the tendency of such 
nearby configurations in phase space to diverge exponentially over time. This divergence is 
parametrized by the Lyapunov exponent for the system, defined as A = \log-g^ where A is 
a distance between two configurations and A is the initial distance at time 0. Numerical 
results shows very fast onset of chaos around the moment t* where the strong turbulence 
begins. 

The highlights of our study for early universe phenomenology are the following. The 
mechanism of preheating after inflation is rather robust and works for many different systems 
of interacting scalars. There is a stage of turbulent classical waves where the initial conditions 
for preheating are erased. Initially, before all the fields have settled into equilibrium with a 
uniform temperature, the reheating temperature may be different in different subgroups of 
fields. The nature of these groupings is determined by the coupling strengths. 
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V. CONCLUSION 

We considered the dynamics of spontaneous symmetry breaking, which occurs when a 
scalar field falls down from the top of its effective potential. We have found [!I| that the 



main part of this process typically completes within a single oscillation of the distribution of 
the scalar field. This is a very unexpected conclusion that may have important cosmological 
implications. 

One of the most efficient mechanisms for the creation of matter after inflation in theories 
with convex effective potentials (V"((j)) > 0) is the mechanism of parametric amplification 
of vacuum fluctuations in the process of homogeneous oscillations of the inflaton field, which 
was called preheating [|l[]. It has also been noted that in the case where potentials become 
concave (V"(</>) < 0), preheating may become more efficient [TJ3|]. Now we see that this effect 



is very generic. In many theories with concave potentials the energy of an unstable vacuum 
state is transferred to the energy of inhomogeneous classical waves of scalar fields within a 
single oscillation of the field distribution. We emphasize here that we are talking about the 
oscillations of the field distribution rather than about the oscillations of a homogeneous field 
4> because quite often the homogeneous component ((f)) of the field remains zero during 
the process of spontaneous symmetry breaking. 

One of the important consequences of our results is the observation |TD[ that in many 
models of hybrid inflation |5j the first stage of reheating occurs not due to homogeneous 
oscillations of the scalar field but due to tachyonic preheating |TD[. In particular, this 



significantly alters the theory of fermionic preheating |11| for the hybrid inflation 



The process of preheating and symmetry breaking may take an especially unusual form 
in the theory of brane inflation |7| based on the hybrid inflation scenario and the mechanism 
of tachyon condensation on the brane antibrane system [^(J . 

The situation in models of the type used in the new inflation scenario is somewhat more 
complicated. In these models the potential is also concave. However, the expansion of the 
universe stretches inhomogeneities of the field rolling down from the top of the effective 
potential and makes it homogeneous on an exponentially large scale. Therefore to evaluate 
a possible significance of tachyonic instability in this regime one must compare the amplitude 
of the homogeneous component of the field with the amplitude of the quantum fluctuations. 
The result appears to be very sensitive to the scale of spontaneous symmetry breaking in such 
models. A preliminary investigation of this issue indicates that in small-field models where 
the scale of spontaneous symmetry breaking is much smaller than M p , the leading mechanism 
of preheating typically is tachyonic. If correct, this would be a very interesting conclusion 
indicating that in large-field models the leading mechanism of preheating typically is related 
to parametric resonance, whereas in small-field models the main mechanism of preheating 
is typically tachyonic, at least at the first stages of the process. 

Finally we should mention that an interesting application of our methods can be found in 
the recently proposed ekpyrotic and pyrotechnic scenario [^1],^]. Even though we are very 
skeptical with respect to the ekpyrotic/pyrotechnic scenario for many reasons explained in 
|22|| , it is still interesting that the methods developed in the theory of tachyonic preheating 



provide us with a very simple theory of the generation of density perturbations in these 
models [^ . 
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